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Abstract. - We present a band structure study of orbital polarization and ordering in the 
two-dimensional triangular lattice transition metal compounds LiV02 and NaTi02 . It is found 
that while in NaTi02 the degeneracy of t2g orbitals is lifted due to the trigonal symmetry of 
the crystal and the strong on cite Coulomb interaction, in LiV02 orbital degeneracy remains 
and orbital ordering corresponding to the trimerization of the two-dimensional lattice develops. 



It is well known that transition metal compounds with orbital degeneracy will in some 
way restructure to remove that orbital degeneracy in the ground state. Well known is the 
example of a two-fold orbitally degenerate case of divalent Cu in octahedral symmetry with 
one hole in a e^-like orbital. A similar case is trivalent Mn in Oh symmetry as in the now 
well known coUossal magnetoresistance materials. In these so called strong Jan Teller systems 
local lattice distortions determine the type of orbital ordering. It is also well established 
that the relative spatial orientation of occupied orbitals on neighboring ions determines not 
only the magnitude but also the sign of the exchange interactions governing the magnetic 
structure of the system In the early 3c? transition metal compounds only the t2g orbitals 
are occupied leaving us with three-fold degeneracy in the cases of Ti'^"'' and V"^^ with one and 
two 3c? electrons respectively assuming also Oh symmetry. In contrast to the Cg orbitals the 
bonding to the neighboring O 2p orbitals is much weaker, bandwidths are much smaller and 
therefore the removal of the orbital degeneracy may be more subtle. It has for example recently 
been suggested that the orbital degeneracy in LiV02 can be lifted by a particular kind of 
orbital ordering driven by the nearest neighbor exchange interactions . The orbital ordering 
proposed there is one which simultaniously removes the frustration in the spin Hamiltonian 
of this triangular two-dimensional lattice and results in a non magnetic singlet ground state. 
Another much discussed two-dimensional triangular lattice spin system is NaTi02 with spin 
i per Ti ion, which in case of uniform Heisenberg interaction would be a frustrated magnetic 
system, making it one of the few remaining possible examples of a resonating valence bond 
(RVB) ground state However also here the orbital degeneracy could have been lifted by 
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Fig. 1. - The fragment of the crystal structure of NaTi02 and LiV02- x,y,z are along the cubic 
crystal axes, the z' axis is the trigonal threefold rotation axis. 



orbital ordering resulting in a non-uniform Heisenberg spin Hamiltonian and a subsequent 
removal of the frustration. In this paper we present a detailed band structure study of these 
two systems both with and without the inclusion of the on site d-d Coulomb interaction in 
a so called LDA+U approximation. We find that although LiV02 behaves as suggested by 
Pen et al ^ NaTi02 should not be considered as orbitally degenerate and remains a candidate 
for an RVB ground state 

LiV02 and NaTi02 crystallize in an ordered rock salt-like structure with alternating [111] 
planes of Li(Na),0,V(Ti),0 ions as shown in fig. |^. Each of these layers forms a triangular two 
dimensional lattice and since the layers containing the magnetic ions are so well seperated the 
exchange interactions between layers is very small compared to that within a layer. Peculiar in 
these compounds is a change from a high temperature paramagnetic phase and a Curie- Weiss 
susceptibility with a large negative effective Curie temperature (which corresponds to strong 
antiferromagnetic coupHng between local moments) to a low temperature nonmagnetic system 
without any sign of long-range order This is usually explained by the frustration of the 

triangular lattice antiferromagnet leading to the possibility of an RVB ground state with a 
quantum liquid of randomly distributed spin singlet pairs |^ . However such models do not take 
into account orbital degeneracy and the possibility of orbital ordering, which can qualitatively 
change the type of exchange interactions as discussed above. 

While the nearest neighbors of the transition metal ion (oxygen atoms) are arranged in an 
almost perfect octahedron with Oh symmetry, the overall symmetry of the crystal structure is 
trigonal, and in that case the t2g level is split into a nondegenerate Aig and double degenerate 
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Fig. 2. - Left panel - the total densities of states from LDA calculations. Right panel - the partial 
densities of states from LDA calculations: dotted line - Aig, solid line - Eg. The Fermi energy is the 
zero energy. 



Eg levels. (Note: in further discussion eg and Eg have different meanings). The t2g-orbitals can 
be represented as the set {dxy,dyz,dzx) in the coordinate system with axes pointing towards 
the 0-neighbors (fig. 0). In this coordinate system the Aig and iJg-orbitals have the following 
form §]: 

Aig : {dxy + dyz + dzx)/V^ Eg-. {1/ ^/i){dzx - dy^) 

{l/V6){-2dxy+dy,+d,x) 

The resulting orbital order will strongly depend on how this trigonal splitting will affect the 
resulting orbital polarization of the partially filled t2g band. If the centre of gravity of the 
Aig band is lower than Eg, then in the case of the d^ configuration the degeneracy will be 
lifted, but the d^ configuration will be still degenerate provided that the Hunds rule exchange 
is larger than this splitting which is certainly expected to be the case. If the Eg level is lower 
than Aig, then the d^ configuration in the high spin state will be nondegenerate. 

In order to look into this problem we have performed band structure calculations for LiV02 
and NaTi02 by LMTO method Q in standard Local Density Approximation (LDA) Q and 
also in LDA+U approach ^ which allows to take into account Coulomb correlations between 
3d electrons of transition metal ions. We used the crystal structures and atomic positions as 
given by The total DOS obtained in LDA and presented in the left panel of fig. ^ 

shows the O 2p band between -8 and -4 eV, (Ti,V) 3d band of t2g symmetry crossing the Fermi 
level between -1 and 1 eV, and of Cg symmetry around 3 eV well above the Fermi energy. As 




Fig. 3. - (a) - the Ti(3d) partial densities of states for rhombohedral NaTi02 obtained from LDA+U 
calculations, (b) - The V(3d) partial densities of states for LiV02 from LDA+U calculations. The 
upper panel corresponds to the V atom with total spin direction "up", the middle one for the V atom 
with "spin-down", and the lower one to the V atom with zero moment. Only xz and yz orbitals in 
local coordinate system for each V atom of the t2g subband are depicted. On every V-atom z axis 
points to the oxygen atom above to the center of the V-triangle. The Fermi energy is the zero energy. 



one can see, indeed the t2g-eg crystal field splitting is larger than the bandwidth and V 3d 
subbands of the t2g and eg symmetry are separated by a gap. We should mention that the eg 
band width is even smaller than that of t2g, this is because the e^-Cg hopping of the neighboring 
tranisition metal ions goes through oxygen atoms and the angle of Me-O-Mc bond is nearly 
90°, so the Gg-Cg hybridization is small. In spite of the fact that there is a trigonal distortion of 
the lattice the t2g band is not splitted due to the trigonal symmetry into Aig and Eg subbands 
and a more delicate quantitative analysis is needed to clarify this problem. 

In the right panel of fig. || the partial DOS for decomposition of the t2g band into orbitals of 
the Aig and Eg symmetry are presented. One can see that the situation can not be described 
in the simple terms of Aig-Eg " splitting" : both curves have the same width and they are 
approximately in the same energy region. We can estimate the actual splitting of the Aig and 
Eg levels by calculating the values of the centres of gravity of these bands. 

In the case of NaTi02 with one d electron the center of gravity of the Aig band is 0.1 eV 
lower than that of the Eg band. As a result the occupied part of t2g band has slightly more Aig 
character than Eg, and the occupancy of orbitals are .25 and .20 per spin-orbital for Aig and 
for each Eg correspondingly. This means that the degeneracy of the t2g orbitals is essentially 
lifted but the splitting is still much less than the band widths. This small splitting is non the 
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Fig. 4. - The t2g holes from LDA+U calculations for LiV02. Only V atoms for one triangle in a 
hexagonal plane are drawn. The view is from the point directly above V-triangle 



less important since as we will see below if we turn on the d-d Coulomb interaction in LDA+U 
the Aig band will be occupied and the Eg unocupied now with a splitting mainly due to U. 
However the choice as to which band is occupied and which one not is dictated by the small 
crystal field splitting. The above result would also indicate that the Aig-Eg local excitation 
energy would be only 0.1 eV or so and would contribute to charge conserving excitonic-like 
excitations. In this LDA+U calculation the d-d Coulomb interaction was found to be 3.6 eV 
(taking into account the screening of t2g electrons by eg electrons ) which is much larger 
than the t2g band width and leads to the localization of a single d electron in the Aig orbital. 

The LDA+U calculations were carried out for both antiferromagnetic and ferromagnetic 
cases. For the AFM case we choose the simplest magnetic order with four nearest neighbors 
out of the six in the basal plane having anti-parallel spin orientation and other two parallel. 
Independent of the spin ordering a single d electron in the t2g shell of the Ti ion turned out 
to be localized in the Aig orbital. The occupation numbers for the majority spin are 0.9 
for Aig and 0.1 for Eg for both FM and AFM cases. The Ti{3d) projected density-of-states 
obtained from LDA+U calculations is shown in fig. ^(a). So we can say from the fig. |^(a) that 
the LDA+U solution for NaTi02 is almost fully orbitally polarized. One can see from eq. |^ 
that the Aig orbital {d^z^^r^ in fig. §(a)) is symmetric in the hexagonal Ti-Ti plane and the 
occupation of this orbital leads to the isotropic exchange. This indicates that NaTi02 would 
still behave like a frustrated spin system. 

In the case of LiV02 the centre of gravity of the Aig band is only 0.025 eV lower then the 
centre of gravity of the Eg band, and the resulting occupancies are 0.37 and 0.36 for Aig and 
for each i?g-orbitals correspondingly. In this situation orbital degeneracy is not lifted, since we 
now have two electrons one in a Aig and one in a Eg orbital, because the Hunds rule exchange 
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strongly favours the high spin state. As a result the appearance of some kind of orbital order 
can be expected. In |^,|l3|,|l^ the formation of local spin singlets on trimers containing V-atom 
triangles was suggested as the model explaining the low-temperature nonmagnetic behavior of 
LiV02 . Those spin singlets were stabilized by a specific orbital order j2| . 

The LDA+U method is based on a mean-field approximation and can not fully reproduce 
the essentially many-electron singlet wave function, especially the correct energy difference of 
singlet-triplet configurations. However a single Slater determinant trial wave function can still 
describe the basic relationship between spin and orbital degrees of freedom. To imitate trimer 
spin singlets we performed LDA+U calculations with spin-order of the type "up-down-zero" 
on every triangle (closed circles on the fig. |l|). The self-consistent calculations resulted in the 
orbital order of the same type as proposed in 1^ from model calculations [fig. |^(b)]: on every 

V atom the occupied orbitals are xz and yz if in a local coordinate system z axis is directed 
towards the oxygen atom sitting just above the center of V-triangle, and x and y axes are 
directed towards other oxygens of an octahedron(fig. |l|). We should also mention the fact that 
the LDA-I-U calculations give the correct semiconducting state for LiV02 [fig. §(b)] instead of 
a metallic state from "normal" LDA (left panel of fig. ||). 

In fig. ^ the angular distribution of the t2g hole is presented as was obtained from the 
LDA+U calculations. It indicates the same orbital order proposed in [^ (fig. 1(a) in ||]): xz 
and yz orbitals are occupied, the t2g hole is in xy orbital in a local coordinate system of every 

V atom. 

Both LiV02 and NaTi02 were regarded as candidates for realization of Anderson's "res- 
onating valence bond" systems with a quantum liquid of randomly distributed spin singlet 
pairs. Our results show that while in LiV02 more complicated trimer spin singlets with 
corresponding orbital order are formed, no orbital order due to the crystal field lifting of 
the orbital degeneracy is present in NaTi02 , and its magnetic properties are most probably 
explained by a nondegenerate model, so that it is indeed a good candidate for Anderson's 
RVB state. What then is the nature of the structural phase transition observed in NaTi02 at 
Tc=250, remains an open question. 

Summarizing, we have shown that the degeneracy of the <2g-orbitals in NaTi02 is lifted 
because of the trigonal symmetry of the crystal and the large d-d Coulomb interaction and 
no orbital ordering occurs. In LiV02 orbital degeneracy remains in spite of the same trigonal 
distortion as in NaTi02 , and in effect the orbital ordering consistent with the trimerization 
of the two-dimensional lattice takes place. 
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